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We clarify the origin of anomalous tunneling [Yu. Kagan et al. Phys. Rev. Lett. 90 (2003) 
130402] i.e. the perfect transmission at low energy limit of tunneling of phonon excitations 
across the potential barrier separating two Bose condensates. The perfect transmission is a 
consequence of the coincidence of the wave function of the excited state at low energy limit 
and the macroscopic wave function of the condensate. We show that the perfect transmission 
at low energy occurs even at finite temperatures within the scheme of Popov approximation. 
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Bose-Einstein condensation in trapped dilute alkali- 
metal atoms has opened up a chance to find new quan- 
tum phenomena, owing to the controllability of the 
strength of interaction by Feshbach resonance and of a 
one-body potential by focused laser beam. With these 
backgrounds, many theoretical studies have been done 
with the Gross-Pitaevskii(GP) and Bogoliubov equa- 
tions,"'^ which are no longer theoretical toy models but 
describe weakly interacting Bose gases realized in alkali- 
metal atoms. One of those theoretical studies has been 
done by Kovrizhin^ and Kagan et al.^ They considered 
tunneling of excitation over the barrier separating the 
two condensates and found that the perfect transmission 
occurs at low energy limit, which was named "anomalous 
tunneling" Subsequently, Danshita et al.'*"'' considered, 
extensively, related problems such as tunneling between 
two condensates with different phases^"® and the excita- 
tion spectrum in Kronig-Penny potential.^'* The trans- 
mission of excitations through randomly distributed po- 
tential barriers was discussed by Bilas and Pavloff.^ 

Kagan et al. attributed the anomalous tunneling to 
a virtual resonance level at low energy, which results 
from the repulsive potential barrier and the depletion of 
condensate wave function near the barrier. Through the 
study of tunneling between the two condensates with dif- 
ferent phases, on the other hand, Danshita et al.^ found 
that anomalous tunneling occurs only when the local- 
ized component is nonzero in the asymptotic solution 
of Bogoliubov equation. Earlier arguments^''' of the ori- 
gin of anomalous tunneling, however, rely on the explicit 
analytical solutions of GP and Bogoliubov equations in 
the presence of a particular form of potential barrier; 
V{x) oc S{x) or the rectangular well V{x) — Vb(> 0) for 
|x| < a and otherwise (x-axis is taken to be the normal 
direction of the wall of the potential barrier) . Therefore 
their arguments can not tell whether anomalous tunnel- 
ing occurs or not in the presence of general form of po- 
tential barrier, which does not allow the exact analytical 
solution for GP and Bogoliubov equations. Furthermore, 
needed is an interpretation applicable to finite tempera- 
tures, where the exact analytical solution is not available 
in general. 
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Fig. 1. Set-up of the system. The two condensates are separated 
by the barrier with the shape of wall. The normal direction is 
taken to be 2:-axis. The incident angle of an excitation to the 
wall is denoted by </>. 



In this Letter, we show the anomalous tunneling occurs 
for the potential barrier V{x) being arbitrary but short- 
ranged and even function (V{—x) — V{x)). Further 
within the Popov approximation, we show the anoma- 
lous tunneling occurs even at finite temperatures. 

First we consider weakly interacting Bose system at 
zero temperature in the presence of potential barrier 
V{x). Linear dimensions in the transverse directions 
(y, z) are assumed to be large so that the mean-field ap- 
proximation is valid. The system we consider is schemat- 
ically described in Fig. 1. We take the incident angle (p to 
be arbitrary within the range [0,7r/2), while has been 
taken to be zero in earlier studies.^"'* 

In the mean-field theory, the macroscopic wave func- 
tion \l/ satisfies the stationary Gross-Pitaevskii equation 



with 



= 0, 



H = — + V{x) - n + C/|*p 



2m 



(1) 



(2) 



Here U is the coupling constant of two-body short-range 
interaction and is related to Airh^aQ/m with the scat- 
tering length oq. Other notations are conventional. The 
system is uniform in j/, z directions and hence ^' depends 
only on x. Far away from the barrier 3> 1, 5* ap- 
proaches y/JiJU . The spatial variation of "^{x) is scaled 
by the healing length ^ = h/ ^Jrnfl. Henceforth "^{x) is 
considered to be real. In the presence of V{x) and ^'(a;), 
the two-component wave function (u(r; e), v{r\ e))* of the 
state with excitation energy e satisfies the Bogoliubov 
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equation 

H' -U^^{x) 
f/**2(a;) -H' 



u{r; e) 
v{r; e) 



u{r; e) 
v{r; e) 



function. When the solution of (7) has the asymptotic 
form: 



with H' = + Owing to the translational 

invariance in y, z directions, the solution of (3) has the 
form of 



(e"=ii^ + r(e)e-"=ii^) 
(3) 

for X <^ —I and 

1 



1 



1 



t{e)e 



ik\\x 



+ c(e)e" 



1 



(13) 



(14) 



u{r; e) 
v{r; e) 



cxp(ifcj,y + ikzz) 



u{x; e) 

v{x; s) 



(4) 



The wave function {u{x;e),v{x;e)Y satisfies the follow- 
ing equation 



H" 



-U<b'^{x) 
-H" 



u{x; s) 
v{x;e) 



u{x; e) 
v{x;e) 



with 



H" = 



and k\ 



2m dx^ 



1,2 I 1,2 



2m 



V{x)- i^ + 2U\'^f 



Intro ducing reduced variables x = x/S,, ^{x) = 
^/U/ fj,^{x), V = V/ji, k±_ = k±^ and s = e/ji, equa- 
tions (1) and (5) reduce to the following equations: 



with/i = -i^ 



m{x) = 

+ V{x)- l + ^'^{x) and 



(6) 



for X ^ I, t{e) and r(e) arc the transmission and re- 
flection coefficients of excitation with energy e for the 
potential barrier V. 

Now wc derive the solution of (7) corresponding to (13) 
and (14) in another way. In terms of 

S{x; e) = u{x; e) + v{x; e) 

and 
(5) 

G{x; e) = u{x; e) — v{x; e), 

(7) becomes 

{h + kl/2)S{x;e) = sG{x;s), (15) 

Ch^-'> + kl/2)G{x;e) = eS{x;e), (16) 

with h^-'> = h + 2*2(a;). Expanding S{x;e) and G{x;e) 
as 

oo oo 

S{x;e) = ^e^S'^^Hx), G{x;s) = ^£"G("Hx) 



n=0 



n=0 



h" -*2(s) 
*2(x) -h" 



v{x; s) 



u{x; s) 
v{x; e) 



(7) 



(17) 

with respect to e, the equations for S^^\x) G^^\x) and 
S^^\x) are obtained as 

hS^°\x) = 0, h^-^G^°\x) = 0, (18) 



with h" = h + ^"^{x) + k\/2, respectively. Henceforth, 
we omit the bar for simplicity. At |a;| ;:?> 1, where F ~ 
and ^{x) ~ 1, the basis of the solution of (7) is given by 
the two plane-wave solutions 



„ibife||X 



and exponentially growing or converging solutions 



1 

6.{e) 



(8) 



(9) 



with e±{e) = ±s/l+^ - £. The wave number A;|| in the 
plane-wave solutions (8) is given by 



and 

hS^^\x) = G^^\x). (19) 

Equations (18) and (19) arc sufficient to derive the low 
energy properties of tunneling of excitations over the po- 
tential barrier. 

The solution of the first equation of (18) can be written 
as the linear combination 

5(°)(,T)=4"),Se(.T)+4")so(x) (20) 

of the even function Se{x) and odd function So{x) satis- 
fying 

dso(x) 



hSe{x) — 0, 



3(0) 



1, 



^11 ~ y ^ ^1' with k = 
It is convenient to express 

A;||=A;cos^, k± = ksm(f>, (/)G[0,7r/2) 



2(^1 + ^2-1). (10) and 



hs,{x) = 0, So(0) = 0, 



da; 

dso(x) 







x=0 



(11) 



dx 



-1, 



(21) 



(22) 



respectively. At |x| ^ 1, where V"(x) ^ and 'i'(x) ^ 1, 



in terms of the incident angle with respect to the wall h reduces to -ig^ and as a result, Se{x) and So{x) have 
of the potential barrier (see Fig. 1). 

The growing rate (or converging rate) k in (9) is given 

by 



the forms of 

Se{x) ae + I3e\x\, So{x) ~ aoSgn(a;) -|- PoX, (23) 



At=y2(7 



l + e^ + l) + k^ sin^ 



In (9), the exponentially-diverging component 



(12) 



should be absent in the physical solution for the wave 



respectively. Generally, etc, cto and /3o arc nonzero; wc 
have confirmed that this is true for V{x) being the rect- 
angular well^ analytically and V{x) = exp(— x^) numer- 
ically. On the other hand, the relation 



/?e = 



(24) 
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holds as a result of 

s,{x) = *(a;)/*(0), 
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(25) 



which follows from the fact that both hand sides of (25) 
satisfy the same equation and initial condition (21) and 
^'(a;)/5'(0) approaches a constant at |a;| » 1. For the 
moment, wc leave the constant f3e in the expressions of 
physical quantities in order to make the role of (24) man- 
ifest in the perfect transmission at low energy limit. 

The solution of the second equation of (18) can be 
written as the linear combination 
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for —1/s <C a; <C — 1 and 

S{x; £) ~ (2 - e)t{s) (1 + iix) + 0{e^) (40) 

for 1 <c a; <C l/e. Here i denotes ecoscj). Equating (35) 
with (39), we obtain 

ie(£)ae - io(£)Qo = 1 + r{e) (41) 



and 



-ie(e)/3c + io(£)/3o = (1 - r(e))ie, (42) 



G(°)(a;) = B(")5e(x) + i?f ffo(a;) 



(26) 



with Ae,o(e) = A,o(£)/(2 - e). Equating (36) with (40), 
we obtain 



of the even function ge{x) and odd function go{x) satis- 
fying 



and 



/i(-)5e(a;) = 0, 5e(0) = l, 



dfl'e(a;) 



and 



h^-^go{x) = 0, 5o(0) = 0, 



da; 
d3o(a;) 



da; 



x=0 



x=0 



(27) 



(28) 



Ae(e)ae + Ao{e)ao = t{e) 



ie(£)/3e + io(e)/3o = yet{e). 



(43) 



(44) 



The four equations (41), (42), (43) and (44) determine 
^e,o(e), ?■(£) and t{e) uniquely. Consequently, we obtain 



t{e) = 



respectively. At |a;| ^ 1, where ^ reduces to — 
fife (a;) and go{x) have the forms of 

g^{x) ~ a^-^e^l^l + /3(-)e-2|^l (29) 

and 



(/3c - ieae){(3o - iecco) 
Setting /3e = (24), we arrive at 



+ O(e'). (45) 



t{E) = ( 1 - i£^ ) + 0(£2) = 1 + i£ cos 0^ + 0{e% 

(46) 



go{x) ~ (a^-^e'l^l + /J^^e-'l^l) sgn(a;) (30) i-^- t^^e Perfect transmission 



respectively. Generally, the coefficients ae,o^ of the expo- 
nentially growing function arc not zero. In (26), on the 
other hand, those exponentially growing terms should be 
absent. It is possible only when = = Q and hence 

G(°)(a;)=0 (31) 

follows. S'^^\x) satisfies 

hS^^\x)=G^'^\x)=Q (32) 

and hence can be written as 

AWse(x)+4i)so(a;). (33) 
With use of it, ^(a;; e) is written as 

S(x; e) = A,{e)s,{x) + A,{e)s,{x) + ©(e^) (34) 
with y4o.o(e) = ^o"o + £^o!"o. S{x; e) becomes 

^c(e)ac-^o(e)ao + (-^(e)/3e+^o(e)/3o)a;+0(e^) (35) 
at .T <§; —1, and 

^e(e)ac + Ao(e)ao + (Ae(£)/3e + ^o(£)/3o)a; + 0(e') (36) 
at a; 1. 

Prom (13), on the other hand, we obtain 
S{x; £) ~ (1 + 6l+(e))(e"=ii^ + r(e)e-'^ii^) 
for a; <C — 1 and 

5(a;;£)~(l + ^+(£))t(e)e''=ii" 



lim tie) = 1 

at low energy limit. From (46), we obtain the phase shift 

5{e) 



(5(e) =e cos 0^+O(£2), 

Po 



(47) 



(37) 



which is defined as {t{e) = \t{e) \ exp(i(5(£))). The energy 
dependence of phase shift (47) is consistent with the ear- 
lier result.'^ 

The anomalous timneling is a direct result of (24), 
which follows from the relation (25) or the fact that the 
even part {x) of the wave function of excitation at low 
energy limit coincides with the wave function of conden- 
sate. For the anomalous tunneling, it is crucial that the 
low energy excitation is obtained by infinitesimal defor- 
mation of the ground state. This applies to the present 
case; the low energy excitations are the Goldstone-mode 
as a consequence of U(l) symmetry breaking. 

With this finding, we discuss the anomalous tunnel- 
ing at finite temperatures. Among several theories for 
Bose system at finite temperatures, we take the Popov 
approximation^ as an example. From now on, we re- 
cover the conventional unit from the dimensionless one. 
In the Popov approximation, the condensate wave func- 
tion ^(a;) satisfies 

/?p*(a;) = (48) 



(38) with 



fora;> 1. With use of ^ e cos(l)+0{e^), k = 2+0{e^) 
and 1 + ^^.(e) = 2 £ + ©(e^), we obtain 

S{x;e) - (2-£) (1 -|-r(£) + {I - r{e))lex)+0{e^) (39) 



2m 



+ V{x) - M + U\^{x)\'' + 2Up{x). 
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The wave function of excitation satisfies 
V nr;e) J 

with H^ = Hp + J7|*(a;)|^. The function p{x) is given by 
p{x) = J2 {Wi{r)\y{s,) + Mr)\\l + /(£,))) . (50) 

i 

Here i is the index of the cigcnstatc for (49) under a 
certain boundary condition and Ui{r) and Vi{r) are the 
wave functions of the i-th eigenstate with energy Si and 
are normaUzed by 

J dr(k(r)|2-k(r)p)=l. 

The function f{e) is given by the Planck distribution 
function {exp[e / {k^T)] — 1)~^ at temperature T. Note 
that in the presence of V{x) and '^{x), \ui{r)\'^ and 
|t;i(r)p in the right-hand side of (50) depend only on 
X. Prom (49), we see that 

S(x) = lim{u(r;e) + v{r;s)) (51) 

£— »0 

satisfies (48) and the even part of (51) is proportional 
to \E'(a;). Following the argument at zero temperature, 
we conclude that anomalous tunneling occurs at finite 
temperatures within the Popov approximation. Indeed, 
we have confirmed numerically that anomalous tunnel- 
ing occurs for V{x) oc exp(— x^) at finite temperatures 
within the Popov approximation.^*^ We expect that the 
anomalous tunneling occurs at finite temperatures even 
when other approximate theories are taken if low energy 
excitations are infinitesimal deformations of the ground 
state. 

So far we have considered the tunneling by scatter- 
ers with the wall geometry. Lastly, we remark on the 
scattering of excitations by scatterers with other geome- 
tries.^^' Applying the same analysis used so far, we 
obtain the cross section a{e) which depends on energy 
as 

cr{e) (X ~ fc^ (52) 



Author Name 

at low energy for the scattering of excitation by a cylin- 
drical potential and 

a(e) cx e'^ ~ k'^, (53) 

for the scattering of excitation by a spherically symmet- 
ric potential. (52) and (53) result from the fact that the 
wave function of excitations in s-wave channel at low en- 
ergy coincides with that of condensate. Note that (53) is 
consistent with the Rayleigh scattering of sound wave in 
classical wave mechanics. ^"^ The method presented in this 
Letter is useful to discuss low energy properties of tun- 
neling and scattering problems by scatterers with various 
geometries in condensed bose systems. 

In conclusion, we have shown that the anomalous tun- 
neling occurs for symmetric potential barrier V{x) = 
V" (—.'/;). The proof for the occurence of the anomalous 
tunneling over asymmetric potential barrier V{x) ^ 
F(— x) is left as a future problem. Another open issue 
is the origin of anomalous tunneling in the presence of a 
condensate current.^ 
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